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We reevaluate the hadronic contributions to the muon magnetic anomaly, and to the running of 
the electromagnetic coupling constant at the Z-boson mass. We include new ir + ir~ cross-section 
data from KLOE, all available multi-hadron data from BABAR, a reestimation of missing low-energy 
contributions using results on cross sections and process dynamics from BABAR, a reevaluation of 
all experimental contributions using the software package HVPTools together with a reanalysis of 
inter-experiment and inter-channel correlations, and a reevaluation of the continuum contributions 
from perturbative QCD at four loops. These improvements lead to a decrease in the hadronic 
contributions with respect to earlier evaluations. For the muon g — 2 we find lowest-order hadronic 
contributions of (692.3±4.2)-10~ 10 and (701.5±4.7) ■ 10 -10 for the e + e"-based and r-based analyses, 
respectively, and full Standard Model predictions that differ by 3.6 a and 2.4 a from the experimental 
value. For the e + e~-based five-quark hadronic contribution to ct(M z ) we find Aa[^ d (M z ) = (275. 7± 
1.0) ■ 10 -4 . The reduced electromagnetic coupling strength at Mz leads to an increase by 7 GeV 
in the central value of the Higgs boson mass obtained by the standard Gfitter fit to electroweak 
precision data. 



I. INTRODUCTION 

The Standard Model (SM) predictions of the anoma- 
lous magnetic moment of the muon, a M , and of the run- 
ning electromagnetic coupling constant, a(s), are limited 
in precision by contributions from virtual hadronic vac- 
uum polarisation. The dominant hadronic terms can be 
calculated with a combination of experimental cross sec- 
tion data, involving e + e~ annihilation to hadrons, and 
perturbative QCD. These are used to evaluate an energy- 
squared dispersion integral, ranging from the ir°j thresh- 
old to infinity. The integration kernels occurring in the 
dispersion relations emphasise low photon virtualities, 
owing to the 1/s descend of the cross section, and, in case 
of a M , to an additional 1/s suppression. In the latter case, 
about 73% of the lowest order hadronic contribution is 
provided by the 7r + 7r~(7) final state, 1 while this chan- 
nel amounts to only 13% of the hadronic contribution to 
a(s) at s = M\. 

In this paper, we reevaluate the lowest-order hadronic 
contribution, a^ ad,LO , to the muon magnetic anomaly, 
and the hadronic contribution, Aahad(-^f ), to the run- 
ning a(Mz) at the Z-boson mass. We include new ir + ir~ 
cross-section data from KLOE [1] and all the available 
multi-hadron data from BABAR [2-9]. We also perform 
a reestimation of missing low-energy contributions us- 
ing results on cross sections and process dynamics from 
BABAR. We reevaluate all the experimental contribu- 
tions using the software package HVPTools [10], includ- 
ing a comprehensive reanalysis of inter-experiment and 



*Now at CERN, CH-1211, Geneva 23, Switzerland. 
1 Throughout this paper, final state photon radiation is implied 
for all hadronic final states. 



inter-channel correlations. Furthermore, we recompute 
the continuum contributions using perturbative QCD at 
four loops [11]. These improvements taken together lead 
to a decrease of the hadronic contributions with respect 
to our earlier evaluation [10], and thus to an accentuation 
of the discrepancy between the SM prediction of a M and 
the experimental result [12]. The reduced electromag- 
netic coupling strength at Mz leads to an increase in the 
most probable value for the Higgs boson mass returned 
by the electroweak fit, thus relaxing the tension with the 
exclusion results from the direct Higgs searches. 



II. NEW INPUT DATA 

The KLOE Collaboration has published new 7r + 7r~7 
cross section data with 7r + 7r~ invariant mass-squared be- 
tween 0.1 and 0.85 GeV 2 [1]. The radiative photon in 
this analysis is required to be detected in the electro- 
magnetic calorimeter, which reduces the selected data 
sample to events with large photon scattering angle (po- 
lar angle between 50° and 130°), and photon energies 
above 20 MeV. The new data are found to be in agree- 
ment with, but less precise than, previously published 
data using small angle photon scattering [13] (supersed- 
ing earlier KLOE data [14]). They hence exhibit the 
known discrepancy, on the p resonance peak and above, 
with other 7r + 7r~ data, in particular those from BABAR, 
obtained using the same ISR technique [2] , and with data 
from t~ — > Tr~n°v T decays [15]. 

Figure 1 shows the available e + e~ — > tt + tt^ cross sec- 
tion measurements in various panels for different centre- 
of-mass energies (y/s). The light shaded (green) band 
indicates the HVPTools average within 1 a errors. The 
deviation between the average and the most precise indi- 
vidual measurements is depicted in Fig. 2. Figure 3 shows 
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FIG. 1: Cross section of e + e — > -k + 7v versus centre-of-mass energy for different energy ranges. Shown are data from TOF [17], 
OLYA [18, 19], CMD [18], CMD2 [20, 21], SND [22] DM1 [23], DM2 [24], KLOE [1, 13], and BABAR [2]. The error bars show 
statistical and systematic errors added in quadrature. The light shaded (green) band indicates the HVPTools average within 
1 a errors. 
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FIG. 2: Comparison between individual e + e — > 7r + 7r cross section measurements from BABAR [2], 
KLOE 10 [1], CMD2 03 [20], CMD2 06 [21], SND [22], and the HVPTools average. The error bars show statistical 
errors added in quadrature. 



0.9 1 
sis [GeV] 

KLOE 08 [13], 
and systematic 




■ BABAR 
■ KLOE08 
KLOE10 - 
SND 

CMD2 03 
'U -— CMD2 06. 
■»'■-■- Other exp 



0.4 0.6 



0.8 



1.2 1.4 1.6 1.8 
\ls [GeV] 



FIG. 3: Relative local averaging weight per experiment versus 
centre-of-mass energy in e + e~ — > -k + tt~. See Figs. 1 and 2 
for references. 



KLOE increases towards the <fi mass, hence outperform- 
ing BABAR above -0.8 GeV. For example, at 0.9 GeV 
KLOE's small photon scattering angle data [13] have sta- 
tistical errors of 0.5%, which is twice smaller than that of 
BABAR (renormalising BABAR to the 2.75 times larger 
KLOE bins at that energy). Conversely, at 0.6 GeV the 
comparison reads 1.2% (KLOE) versus 0.5% (BABAR, 
again given in KLOE bins which are about 4.2 times 
larger than for BABAR at that energy). The discrepancy 
between the BABAR and KLOE data sets above 0.8 GeV 
causes error rescaling in their average, and hence loss of 
precision. The group of experiments labelled "other exp" 
in Fig. 3 corresponds to older data with incomplete ra- 
diative corrections. Their weights are small throughout 
the entire energy domain. The computation of the dis- 
persion integral over the full tt + tt~ spectrum requires to 
extend the available data to the region between thresh- 
old and 0.3 GeV, for which we use a fit as described in 
Ref. [101. 



the weights versus y/ s the different experiments obtain in 
the locally performed average. BABAR and KLOE dom- 
inate the average over the entire energy range. Owing to 
the sharp radiator function, the available statistics for 



We have modified in this work the treatment of the 
w(782) and </>(1020) resonances, using non-resonant data 
from BABAR [3]. While in our earlier analyses, the reso- 
nances were fitted, analytically integrated, and the non- 
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FIG. 4: Cross section of e + e~ — > tt + it~tt versus centre- 
of-mass energy. Shown are data from ND [25], DM1 [26], 
SND [27] CMD [28], CMD2 [29] and BABAR [3]. The error 
bars show statistical and systematic errors added in quadra- 
ture. The light shaded (green) band indicates the HVPTools 
average within 1 a errors. 



tion measurements and HVPTools averages for the chan- 
nels K + K~tt + 7t~ (upper left), 2ir + 2ir~ir° (upper right), 
37T+371-- (lower left), and 27r + 27r"27r° (lower right). The 
BABAR data supercede much less precise measurements 
from M3N, DM1 and DM2. In several occurrences, these 
older measurements overestimate the cross sections in 
comparison with BABAR, which contributes to the re- 
duction in the present evaluation of the hadronic loop 
effects. 

Finally, Fig. 7 shows the charm resonance region above 
the opening of the DD channel. Good agreement be- 
tween the measurements is observed within the given 
errors. While Crystal Ball [47] and BES [48] pub- 
lished bare inclusive cross section results, PLUTO ap- 
plied only radiative corrections [49] following the formal- 
ism of Ref. [50] , which does not include hadronic vacuum 
polarisation. As in previous cases [30] for the treatment 
of missing radiative corrections in older data, we have ap- 
plied this correction and assigned a 50% systematic error 
to it. 



resonant contributions added separately, we now deter- 
mine all the dominant contributions directly from the 
corresponding measurements. Hence the u) and <f> con- 
tributions are included in the ir + ir~ir°, 7r°7, 777, K + K~ , 
K®K° spectra. Small remaining decay modes are consid- 
ered separately. As an example for this procedure, the 
e + e~ — > 7T + 7r~7r cross section measurements, featuring 
dominantly the u) and <f> resonances, are shown in Fig. 4, 
together with the HVPTools average. 

We also include new, preliminary, e + e~ — > 7r + 7r~27r° 
cross section measurements from BABAR [5], which sig- 
nificantly help to constrain a contribution with disparate 
experimental information. The available four-pion mea- 
surements and HVPTools averages are depicted in Fig. 5 
in linear (top) and logarithmic (bottom) ordinates. 2 

The precise BABAR data [6-9] available for several 
higher multiplicity modes with and without kaons (which 
greatly benefit from the excellent particle identification 
capabilities of the BABAR detector) help to discrimi- 
nate between older, less precise and sometimes contra- 
dicting measurements. Figure 6 shows the cross sec- 



2 The new measurements also improve the conserved vector cur- 
rent (CVC) predictions for the corresponding r decays with four 
pions in the final state. Coarse isospin-breaking corrections 
with 100% uncertainty are applied [16]. We find £?cvc( T ~ — ► 
7r-37r°t/ T ) = (1.07 ± 0.06)%, to be compared to the world 
average of the direct measurements (1.04 ± 0.07)% [51], and 
#CVc(t~ -> 2tt-tt + tt°v t ) = (3.79 ± 0.21)%, to be compared to 
the direct measurement (4.48 ± 0.06)%. The deviation between 
prediction and measurement in the latter channel amounts to 
3.2 <r, compared to 3.6 <r without the BABAR data [30]. It is due 
to a discrepancy in mainly the normalisation of the correspond- 
ing r and e + e — spectral functions. It is therefore important that 
the BABAR and Belle experiments also perform these r branch- 
ing fraction measurements as independent cross checks. 



III. MISSING HADRONIC CHANNELS 

Several five and six-pion modes involving 7r°'s, as well 
as KK[nir] final states are still unmeasured. Owing to 
isospin invariance, their contributions can be related to 
those of known channels. The new BABAR cross section 
data and results on process dynamics thereby allow more 
stringent constraints of the unknown contributions than 
the ones obtained in our previous analyses [30, 31]. 

Pais Isospin Classes. Pais introduced [52] a clas- 
sification of A^-pion states with total isospin I = 0,1. 
The basis of isospin wave functions of a given state be- 
longs to irreducible representations of the corresponding 
symmetry group, which are characterised by three integer 
quantum numbers {partitions of N) N\, N 2 , N3, obeying 
the relations Aq > N 2 > N 3 > and Aq + N 2 + N 3 = N. 
The total isospin / is determined uniquely to be / = 
if Aq — N 2 and N 2 — N3 are both even, and / = 1 in 
the other cases. States {Ni, N 2 , N 3 } are composed by 
N 3 isoscalar three-pion subsystems, N 2 — N 3 isovector 
two-pion subsystems, and Aq — N 2 isovector single pions. 

Simple examples are {110} for e + e~ — > tt + it~ (I = 1, 
/9-like), and {111} for e+e~ 7r + 7r~7r° (I = 0, w-like). 
For four pions there are two channels and two isospin 
classes, related at the cross section level by 3 

a{e+e- ^2n + 2n-) = ^310, (1) 

5 

cr(e+e~ -> 7r + 7T~27r°) = -o- 310 + <r 211 . (2) 



3 In the following, and if not otherwise stated, c(X) denotes 
cr(e+e- -¥ X). 
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FIG. 5: Cross section versus centre-of-mass energy of e + e~ — > 2-k + 2ty~ (left) and e + e~ — > tt + ty~2tt° (right), and for linear 
(top) and logarithmic ordinates (bottom). The open circles show data from BABAR [4, 5], which dominate in precision. The 
other measurements shown are taken for the four charged pions final state from ND [32], M3N [33], MEA [34], CMD [35], 
DM1 [36, 37], DM2 [38-40], OLYA [41], CMD2 [42] and SND [43], and for the mixed charged and neutral state from ND [32], 
M3N [44], DM2 [38-40], OLYA [45], and SND [43]. The error bars show statistical and systematic errors added in quadrature. 
The shaded (green) bands give the HVPTools averages. 



The two isospin classes correspond to the resonant final 
states pTTir for {310} and um° for {211}. 

The 1=1 states produced in e + e~ are related to 
the vector part of specific r decays by isospin symme- 
try (CVC). 

Five-Pion Channels. There are two five-pion fi- 
nal states, 27T + 27r~7r° and Tr + ir~3ir°, of which only the 
first has been measured. There is only one isospin class 
{311}, corresponding to utttt and obeying the relation 
cr(27r+27r-7r°) = 2cr(7r+7r-37r°) = §<7 3U . 

BABAR has shown [6] that the first channel is indeed 
dominated by ujtttt, with some contribution from rjinr 
via the isospin- violating decay rj — > 7r + 7r~7r°. These rj 
contributions must be subtracted and treated separately 
as they do not obey the Pais classification rules. At larger 
masses there is evidence for a pirn component, which 
should correspond to 1^ = contributions above the uj. 
Isospin symmetry holds for this contribution. 

The estimation procedure for the unknown five-pion 
contribution is as follows: cr(27r + 27r _ 7r°)„_ exc i = 



cr(27r + 27r 7r°) — a^TT+n ) x B[r\ — > n+n 7T°), 
with B(j] n+n-TT ) = 0.2274 ± 0.0028 [51], 

c^TT+Tr-STr^exc! = ±a(2ir + 2ir-ir°) v _ cxch and 
<7(777r + 7r~) is considered separately. There is no 
contribution from a(rj2Tr°) 1 and the contribu- 
tion of com: with non purely pionic uj decays is 
taken from |er(u;7r + 7r~) x S(w-non-pionic) with 
B(w-non-pionic) = 0.093 ± 0.007 [51]. 

Six-Pion Channels. There are three channels and 
four isospin classes with the relations (1 = 1): 

24 3 

cr(37r + 37r~) = — <7 510 + -CT330 , (3) 

8 2 2 

cr(2-7r + 27T - 27r ) = ^rCsio + g^n + -O330 + 0321 , (4) 

3 3 

O-(7r + 7T~47T ) = 7^510 + -CT411 , (5) 

where the lowest-mass resonant states are p4ir for {510}, 
w3?r for {411}, 3,0 for {330}, and ujpir for {321}. 
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FIG. 6: Cross section data for the final states: K + K~ir + ir~ (upper left), 2-k + 2ty~ tt° (upper right), 3tv + 3tv~ (lower left) and 
27r + 27r _ 27r° (lower right). The BABAR data points are taken from Refs. [4, 6, 7]. All the other measurements are referenced 
in [30, 31]. The shaded (green) bands give the HVPTools averages within la errors, locally rescaled in case of incompatibilities. 



BABAR has measured [7] the cross sections (3) and 
(4), and observed only one p state per event in the 
fully charged mode, thus favouring {510} over {330} 
in (3). The e + e _ — > 2ir + '2ir~ 2tt° process is dominated 
by u!ir + ir~ir° up to 2 GeV. A small ij contribution is also 
found, but only the cross section for -quj is given. 

To estimate the cross section (5) the relative contribu- 
tions of {321} and {411} need to be known, which can 
be constrained from r data. The corresponding isospin 
relations for the r spectral functions are 

V(T~ -> 27T+ 371-^71-°^) = ^"510 + gf411 (6) 

1 1 

+ -=V330 + -V 32 l , 

5 2 

v(t~ -> 7r + 27r~37r ( V) = — u 5 i + -v 4 ii (7) 

4 1 

+ c ^330 + 77^321 , 

5 2 



v(t — > 7T 5tT U T ) 



35 



V51Q ■ 



(8) 



The branching fractions of the first two modes have been 
measured by CLEO. As for the e + e~ final states, they 



are dominated by ui production, oji{ + 2it and coir 2ir°, 
with branching fractions (1.20 ± 0.22) • 10~ 4 and (1.4 ± 
0.5) • 10 -4 , respectively, to be compared to total branch- 
ing fractions of (1.40±0.29) • 10~ 4 and (1.83±0.50) • 10~ 4 
(77 subtracted). This yields the bound U411/U321 < 0.42. 
The limit for ir + ir~4ir° is looser than that quoted in 
Ref. [30], where the {411} partition was assumed to be 
negligible. 

The estimation procedure for the missing six- 
pion mode is as follows: rT(27r + 27r _ 27r )„_ exc i = 
cr(27r+27r-27T ) - <t(t]uj) x B(rj -> n+ir-ir°) x B(ui -> 
ir+n-TT ), with B(u -> ir+n-ir ) = 0.892 ± 
0.007 [51], and a(n + n- 4vr°) = 0.0625cr(37r+37r-) + 
0.145cr(27r+27r-27r°) r) _ excl ± 100%; a (tjit+tt-tt ) = 
a(rjuj) x B{(jJ — > 7r + 7r~7r°) is treated separately, 
and the contribution from non-pionic u> decays 
is given by (1.145 ± 0.145) x a {2^+2^-2^) x 
S(w-non-pionic)/S(cj — > 7r + 7r~7r°). 

KKtv Channels. The measured final states are 
RV^ttT and K+R-jt , with K° s K°n° missing (C Kg = 
— 1). Except for a very small (jm° contribution, these 
processes are governed by K°K*° (890) (dominant) and 
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FIG. 7: Inclusive hadronic cross section versus centre-of-mass 
energy above the DD threshold. The measurements are taken 
from PLUTO [46], Crystal Ball [47] and BES [48]. The light 
shaded (green) band indicates the HVPTools average within 
1 a errors. 



#±#^(890) transitions below 2 GeV. Both I = 0, 1 am- 
plitudes (^4o,i) contribute. The fit of the Dalitz plot in 
the first channel yields the moduli of the two amplitudes 
and their relative phase as a function of mass. Hence 
everything is determined, as seen from the following re- 
lations (labels written in the order KK* with the given 
K* decay modes): 



{K+K-ir° 4- 


-K-R+ir ) = 


> 


-A^ 2 , 


(9) 




+ k2*2«°) = 


> 


+ A,\ 2 , 


(10) 


(K°K-tt+ -1 


-K°K + tt-) = 




+ A,\\ 


(11) 


(K+K°7r- 4 


-K-K°tt + ) = 


\\A 


-Ail 2 . 


(12) 



The measured K^K ± n T cross section (no ordering here) 
is therefore equal to |[|A | 2 + |^i| 2 ] = |(o"rj + 01), and 
a(KKn) = SaiK^K^ir^) for the dominant KK* part. 
Note that, unlike it was assumed in Ref. [30, 31], in gen- 
eral a(K°K° L Tr°) is not equal to a(K+K--K°). 

The complete KKix contribution is obtained from 
a(KKn) =_3a(K ^K^) + cr(07r°) x B((f> -> KK), with 
B(<f> -> KK) = 0.831 ± 0.003, where contributions from 
non-hadronic <f> decays are neglected, whereas decays to 
7r + 7T~7r are already counted in the multi-pion channels. 

KK2tv Channels. The channels measured by 
BABAR are K+K-jt+tt- and K+K^n [9]. They are 
dominated by K*Kir, with Kir not in a K*, and smaller 
contributions from K + K~ p° and (pww. 

In the dominant K*Kw mode one can have 1 = and 
1=1 amplitudes. The different charge configurations 
can be obtained via Ik-k = 1/2 and 3/2 amplitudes, 
where, however, I Klx = 3/2 is not favoured because it 
would have predicted <r(K + K~7r + iT~) = <r(K + K~2ir°), 



whereas a ratio of roughly 4:1 has been measured [9]. 
In the following we assume a pure Ik-k = 1/2 state, so 
that the relevant cross sections read (labels in the order 
K*Ktt, appropriately summing over K°(K )) 

1 



aiK^K^n 
a{K°TT ± K^'K 



a(K°Tr°K 7r 



aiK^K -* 

CT(AT ± 7r =F J ft: ± 7r =F 
a(K ± w a K°'n* 
a(K°-K±K Tt* 



18 
1 

9 
1 

18 
1 

9 
1 

9 
2 

9 
1 

9 
2 

9 



Ao-A.W 
|Ao + ^i| 2 
Ao+A^ 2 , 
Ao+A^ 2 , 
Ao+A^ 2 , 
A a -A 1 \ 2 , 
Ao-A.l 2 . 



(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



This leads to a(KKmr) = 9<j(K + K-Tr°TT ) + 
la(K+K-TT+7r-). 

The inclusive a(KKp) cross section is thus obtained as 
follows: get a^n+n-) = 2cr(027r°) and a(K+K-p°) = 
a(K+K-TT+TT-) - a(if* A: ± 7r = F) - o-^tt+tt") x B{<j> -> 
K + K~) (note that the published BABAR cross section 
table for K*°K ± tt t already includes the branching frac- 
tion for K*° —> K ±r K T ). In lack of more information, 
we assume a(KKp) = Aa(K + K~ p°), with a 100% error, 
and obtain o{KKw) = 9[a{K + AT^tt ) - ct(027t )] + 
laiK^K^) + |a(07r+7r-) + 4a{K+ K~ p°). 

KK3iv Channels. BABAR has only measured the 
final state K + K~Tr + ir~ir° [6], which is dominated by 
K + K~u) up to 2 GeV. The channel 4>rj has been mea- 
sured, and the remaining 07r + 7r _ 7r° amplitude is negli- 
gible. The uj dominance does not apply to the missing 
channels K° 'K* vrF tt+ and K^K^ir^ir , but their dy- 
namics (for instance K*) should be seen in the measured 
K + K~ ir + ir~ ir° mode, so it may be small, at least below 
2 GeV. 

The missing channels are estimated as follows: 
aiK+K-Tr+TT-Tr )^^ = a(K + K-tt+tt-tt ) - a^rj) x 
B((f> — > K + K~) x B(n — > 7r + 7r~7r°). We assume, within a 
systematic error of 50%, a(K° l -?r 7r + 7r~7r )„_ exc i = 
a(K + K~-K + n~iT ) rj _ e - KC \, treat cr((f>n) separately, 
and compute the non-pionic oj contribution by 
2<j(K + K~ir + n~ir°) rj _ exc \ x S(w-non-pionic)/S(w — > 
7T + 7r _ 7r°). Contributions from K° K ± ir T ir + tt~ and 
i^° J ft: ± 7r = F27r below 2 GeV are neglected. 



t]4tv Channels. BABAR has measured 
a(i]2ir + 2Tr~) [6], where the 4-7T state has C = — 1, I = 1. 
Because <7(27r + 27r~) rj cr(7r + 7r~27r°), we assume the 
same ratio for the ?y47r process with the same 47r quantum 
numbers. We thus estimate a(i]4ir) = 2<j(i]2ir + 2ir~), 
and assign a systematic error of 25% to it. 
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IV. DATA AVERAGING AND INTEGRATION 

In this work, we have extended the use of HVPTools 4 
to all experimental cross section data used in the com- 
pilation. 5 The main difference of HVPTools with re- 
spect to our earlier software is that it replaces linear in- 
terpolation between adjacent data points ("trapezoidal 
rule") by quadratic interpolation, which is found from 
toy-model analyses, with known truth integrals, to be 
more accurate. The interpolation functions are locally 
averaged between experiments, whereby correlations be- 
tween measurement points of the same experiment and 
among different experiments due to common systematic 
errors are fully taken into account. Incompatible mea- 
surements lead to error rescaling in the local averages, 
using the PDG prescription [51]. 

The errors in the average and in the integration for 
each channel are obtained from large samples of pseudo 
Monte Carlo experiments, by fluctuating all data points 
within errors and along their correlations. The integrals 
of the exclusive channels are then summed up, and the 
error of the sum is obtained by adding quadratically (lin- 
early) all uncorrelated (correlated) errors. 

Common sources of systematic errors also occur be- 
tween measurements of different final state channels and 
must be taken into account when summing up the exclu- 
sive contributions. Such correlations mostly arise from 
luminosity uncertainties, if the data stem from the same 
experimental facility, and from radiative corrections. In 
total eight categories of correlated systematic uncertain- 
ties are distinguished. Among those the most significant 
belong to radiative corrections, which are the same for 
CMD2 and SND, as well as to luminosity determinations 
by BABAR, CMD2 and SND (correlated per experiment 
for different channels, but independent between different 
experiments) . 



V. RESULTS 

A compilation of all contributions to aJj ad:LO and to 
^Q ; had(-^f ), as well as the total results, are given in Ta- 
ble II. The experimental errors are separated into statisti- 
cal, channel-specific systematic, and common systematic 
contributions that are correlated with at least one other 
channel. 

Table I quotes the specific contributions of the var- 
ious e + e~ — > ir + n~ cross section measurements to 
a had,LO_ Also given are the corresponding CVC-based 



TABLE I: Contributions to o£ ad ' LO (middle column) from the 
individual ir + ir~ cross section measurements by BABAR [2], 
KLOE [1, 13], CMD2 [20, 21], and SND [22]. Also given 
are the corresponding CVC predictions of the r — > ir~ir°v T 
branching fraction (right column), corrected for isospin- 
breaking effects [15]. Here the first error is experimental and 
the second estimates the uncertainty in the isospin-breaking 
corrections. The predictions are to be compared with the 
world average of the direct branching fraction measurements 
(25.51 ± 0.09)% [51]. For each experiment, all available data 
in the energy range from threshold to 1.8 GeV (m T for Bcvc) 
are used, and the missing part is completed by the combined 
e + e~ data. The corresponding (integrand dependent) frac- 
tions of the full integrals provided by a given experiment are 
given in parentheses. 



Experiment a„ 



[io- 



BABAR 
KLOE 
CMD2 
SND 



514.1 ±3.8 (1.00) 
503.1 ±7.1 (0.97) 
506.6 ±3.9 (0.89) 
505.1 ±6.7 (0.94) 



Bcvc [%] 

25.15 ±0.18 ±0.22 (1.00) 
24.56 ±0.26 ±0.22 (0.92) 
24.96 ±0.21 ±0.22 (0.96) 
24.82 ±0.30 ±0.22 (0.91) 



r — > ir~ir°v T branching fraction predictions. The largest 
(smallest) discrepancy of 2.7 a (1.2 a) between predic- 
tion and direct measurement is exhibited by KLOE 
(BABAR). It is interesting to note that the four 
a hadxoj 7r + 7r -j determinations in Table I agree within er- 
rors (the overall \ 2 °f their average amounts to 3.2 for 
3 degrees of freedom), whereas significant discrepancies 
are observed in the corresponding spectral functions [10]. 
Since we cannot think of good reasons why systematic 
effects affecting the spectral functions should necessarily 
cancel in the integrals, we refrain from averaging the four 
values with a resulting smaller error. The combined con- 
tribution is instead computed from local averages of the 
spectral function data that are subjected to local error 
rescaling in case of incompatibilities. 

The contributions of the J/ip and tp(2S) resonances 
in Table II are obtained by numerically integrating 
the corresponding undressed 6 Breit-Wigner lineshapes. 
Using instead the narrow-width approximation, <jr — 
12n 2 T° e /M R ■ S(s ~ M|), gives compatible results. The 
errors in the integrals are dominated by the knowledge 
of the corresponding bare electronic width Y Q R ^ eR . 

Sufficiently far from the quark thresholds we use four- 
loop [11] perturbative QCD, including 0(ev|) quark mass 
corrections [53], to compute the inclusive hadronic cross 
section. Non-perturbative contributions at 1.8 GeV were 
determined from data [54] and found to be small. The 
errors of the -Rqcd contributions given in Table II ac- 



4 See Ref. [10] for a more detailed description of the averaging and 
integration procedure developed for HVPTools. 

5 So far [10], only the two-pion and four-pion channels were fully 
evaluated using HVPTools, while all other contributions were 
taken from our previous publications, using less sophisticated 
averaging software. 



The undressing uses the BABAR programme Afkvac correcting 
for both leptonic and hadronic vacuum polarisation effects. The 
correction factors amount to (1 — U(s)) 2 = 0.956 and 0.957 for 
the J/ip and %1>(2S), respectively. 
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TABLE II: Compilation of the exclusive and inclusive contributions to a]} ad ' LO and Acvhad(Aff). Where three (or more) errors 
are given, the first is statistical, the second channel-specific systematic, and the third common systematic, which is correlated 
with at least one other channel. For the contributions computed from QCD, only total errors are given, which include effects 
from the a s uncertainty, the truncation of the perturbative series at four loops, the FOPT vs. CIPT ambiguity (see text), 
and quark mass uncertainties. Apart from the latter uncertainty, all other errors are taken to be fully correlated among the 
various energy regions where QCD is used. The errors in the Breit-Wigner integrals of the narrow resonances J/tf and ip(2S) 
are dominated by the uncertainties in their respective electronic width measurements [51]. The error on the sum (last line) 
is obtained by quadratically adding all statistical and channel-specific systematic errors, and by linearly adding correlated 
systematic errors where applies (see text for details on the treatment of correlations between different channels). 



Channel a^ d - LO [1(T 10 ] A« had (Ml) [1(T 4 ] 



TT°-y 

7T 7 


^i.^i^j zn u.uo zn u.lo zn u.l^ 


IP 4 01 4 01 4 01 
u.ou zn u.u± zn u.U-L zn u.U-L 


Vl 


n R4 4- n n9 4 ni 4 n ni 

u.u^ zn u.u^ zn u.ul zn u.U-L 


OS 4 00 4 00 4 00 
u.uo zn u.uu zn u.uu zn u.uu 




^07 80 + 1 99 4- 9 ^fl + ^ifi 
uu 1 .ou zn l.__ zn zj.uu zn u. ju 


QA A^i 4 07 4 1 7 4 04 
O'l.'lo zn u.ui zn u.i l zn u.U'L 




dp, 00 4 49 4 1 OR + (1 QS 

zn U.^fc^ Z1Z l.UO zn U.J70 


4 ■iS 4 04 4 1 1 4(1 (1Q 
i.uo zn u.u^ zn u.J-j. zn u.uy 


2h + 2it~ 


1 3 4 1 4 4^ 4 9Q 
J-O.oo zn u.-LU zn u.^to zn u.^i? 


idQ+nn^ + n 1 9 4 os 

o.^ty zn u.uo zn u.J-^ zn u.uo 


TX + TX~2n 


18 01401441 1740 40 
lo.ui zn u . -L^± zn l . ± i zn u . ^tu 


A A'i 4 0*3 4 9Q 4 10 
^±.4:0 zn u.uo zn u.^y zn u. iu 


9ir~'~ 9tt ~~ ir° in pypI "1 

_/( £i l\ /I U/ CALL 1 


n 79 4. n f]A 4 07 4 0^ 
u. / zn u.U'L zn u.u ( zn u.uo 


99 4 01 4 09 4 01 
u . _ zn u.U-L zn u.u^ zn u.U-L 


11 /l O/l \ll CALL, LL LJLiL lOUOJJILl 1 


36 4 09 4 03 4 01 

U.OU —1— U .U Li _L_ U.UU _L_ U.U-L 


11 4 01 4 01 4 00 

U . -L -L • U.U-L 1 U.U-L 1 U.UU 


3tv + 3tt~ 


19-1-001 _i_ n ni _i_ n nn 
u.J-^r zn u.uj. ziz u.ui zlz u.uu 


n ha _i_ n nn i n nn in nn 

U.U4: ZLZ U.UU ZLZ U.UU ZLZ U.UU 


97r + 97r — 9ir° In pttcI "1 

£1 l\ £t IK [II CALL 1 


n 7Q _i_ n nc _i_ n c\a _|_ n HQ 
u. 1 u zlz u.uu zlz u.u* 1 * zn u.uy 


n oc _i_ n no _i_ n no _i_ n no 

U.ZiU ZLZ U.Uj£ ZLZ U.U^ ZLZ U.UO 


7r + 7r _ 47r° (rj excl., from isospin) 


0.11 ±0.01 ±0.11 ±0.00 


0.04 ±0.00 ±0.04 ±0.00 




1.15 ±0.06 ±0.08 ±0.03 


0.33 ±0.02 ±0.02 ±0.01 


r\u) 


0.47 ±0.04 ±0.00 ±0.05 


0.15 ±0.01 ±0.00 ±0.02 


r/2n + 2TV~ 


0.02 ±0.01 ±0.00 ±0.00 


0.01 ± 0.00 ± 0.00 ± 0.00 


rj7v + ir~ 2tt° (estimated) 


0.02 ±0.01 ±0.01 ±0.00 


0.01 ±0.00 ±0.00 ±0.00 


luty (uj —> n°j) 


0.89 ±0.02 ±0.06 ±0.02 


0.18 ±0.00 ±0.02 ±0.00 


cuir + ir~ , cj27t° (oj — > ir°y) 


0.08 ±0.00 ±0.01 ±0.00 


0.03 ±0.00 ±0.00 ±0.00 


lu (non-37r, ivy, 777) 


0.36 ±0.00 ±0.01 ±0.00 


0.03 ±0.00 ±0.00 ±0.00 


K+K- 


21.63 ±0.27 ±0.58 ±0.36 


3.13 ±0.04 ±0.08 ±0.05 


K° S K° L 


12.96 ±0.18 ±0.25 ±0.24 


1.75 ±0.02 ±0.03 ±0.03 


4> (non-KK, 3n, 717, 777) 


0.05 ±0.00 ±0.00 ±0.00 


0.01 ±0.00 ±0.00 ±0.00 


KK-k (partly from isospin) 


2.39 ±0.07 ±0.12 ±0.08 


0.76 ± 0.02 ± 0.04 ± 0.02 


KK2n (partly from isospin) 


1.35 ±0.09 ±0.38 ±0.03 


0.48 ±0.03 ±0.14 ±0.01 


KK3n (partly from isospin) 


-0.03 ±0.01 ±0.02 ±0.00 


-0.01 ±0.00 ±0.01 ±0.00 




0.36 ±0.02 ±0.02 ±0.01 


0.13 ±0.01 ±0.01 ±0.00 


uKK (cj -4 7T°7) 


0.00 ±0.00 ±0.00 ±0.00 


0.00 ±0.00 ±0.00 ±0.00 


J/ip (Breit-Wigner integral) 


6.22 ±0.16 


7.03 ±0.18 


tp{2S) (Breit-Wigner integral) 


1.57 ±0.03 


2.50 ±0.04 


.Rdata [3.7 - 5.0 GeV] 


7.29 ±0.05 ±0.30 ±0.00 


15.79 ±0.12 ±0.66 ±0.00 


Rqcd [1.8-3.7 GeV]„ ds 


33.45 ±0.28 


24.27 ±0.19 


Rqcv [5.0 - 9.3 GeV]udsc 


6.86 ±0.04 


34.89 ±0.18 


Rqcd [9.3 - 12.0 GeV] udscb 


1.21 ± 0.01 


15.56 ±0.04 


Rqcd [12.0 - 40.0 GeV] udscb 


1.64 ± 0.01 


77.94 ±0.12 


Rqcd [> 40.0 GeV] udscb 


0.16 ±0.00 


42.70 ±0.06 


Rqcv [> 40.0 GeV]* 


0.00 ±0.00 


-0.72 ±0.01 



Sum 692.3 ± 1.4 ± 3.1 ± 2.4 ± 0.2^, ± 0.3 QC d 274.97 ± 0.17 ± 0.78 ± 0.37 ± 0.18^, ± 0.52 QC d 
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FIG. 8: Inclusive hadronic cross section ratio versus centre-of- 
mass energy in the continuum region below the DD threshold. 
Shown are bare BES data points [48] , with statistical and sys- 
tematic errors added in quadrature, the data average (shaded 
band), and the prediction from massive perturbative QCD 
(solid line — see text). 



count for the uncertainty in a s (we use a g (Af§) = 
0.1193±0.0028 from the fit to the Z hadronic width [55]), 
the truncation of the perturbative series (we use the full 
four-loop contribution as systematic error), the full dif- 
ference between fixed-order perturbation theory (FOPT) 
and, so-called, contour-improved perturbation theory 
(CIPT) [56], as well as quark mass uncertainties (we use 
the values and errors from Ref. [51]). The former three 
errors are taken to be fully correlated between the vari- 
ous energy regions (see Table II), whereas the (smaller) 
quark-mass uncertainties are taken to be uncorrelated. 
Figure 8 shows the comparison between BES data [48] 
and the QCD prediction below the DD threshold be- 
tween 2 and 3.7 GeV. Agreement within errors is found. 7 



Muon magnetic anomaly. Adding all lowest- 
order hadronic contributions together yields the estimate 
(this and all following numbers in this and the next para- 
graph are in units of 10" 



-lrn 



had.LO = 6g2 3 ± 4 2 

h 1 



(21) 



7 To study the transition region between the sum of exclusive mea- 
surements and QCD, we have computed aJj ad ' LO in two narrow 
energy intervals around 1.8 GeV. For the energy interval 1.75— 
1.8 GeV we find (in units of 10" 10 ) 2.74± 0.06 ± 0.21 (statistical 
and systematic errors) for the sum of the exclusive data, and 
2.53±0.03 for perturbative QCD (see text for the contributions to 
the error). For the interval 1.8-2.0 GeV we find 8.28±0.11±0.74 
and 8.31 ± 0.09 for data and QCD, respectively. The excellent 
agreement represents another support for the use of QCD beyond 
1.8 GeV centre-of-mass energy. Comparing the a^ ad,LO predic- 
tions in the energy interval 2—3.7 GeV, we find 26.5 ± 0.2 ± 1.7 
for BES data, and 25.2 ± 0.2 for perturbative QCD. 



HMNT 07 (e + e -based) 

-285 + 51 I — 

JN 09 (eV) 

-299 ±65 I 

Davier et al. 09/1 (t-based) 
-157 + 52 

Davier et al. 09/1 (eV) 

-312 + 51 I • 

Davier et al. 09/2 (eV w/ BABAR) 
-255 ±49 

HLMNT 10 (e + e"w/ BABAR) 

-259 ± 48 I 

DHMZ 10 (x newest) 
-195 ±54 

DHMZ 10 (eV newest) 

-287 ±49 I— 



BNL-E821 (world average) 

0±63 



-700 -600 -500 -400 



-300 -200 



a exp 



-100 



x10" 



FIG. 9: Compilation of recent results for a^ M (in units of 
10" 11 ), subtracted by the central value of the experimen- 
tal average [12, 57]. The shaded vertical band indicates 
the experimental error. The SM predictions are taken from: 
this work (DHMZ 10), HLMNT [58] (e+e - based, including 
BABAR and KLOE 2010 tt+tt" data), Davier et al. 09/1 [15] 
(r-based), Davier et al. 09/1 [15] (e+e~-based, not including 
BABAR tt+tt" data), Davier et al. 09/2 [10] (e + e"-based in- 
cluding BABAR tt+tt" data), HMNT 07 [59] and JN 09 [60] 
(not including BABAR tt+tt" data). 



which is dominated by experimental systematic uncer- 
tainties (cf. Table II for a separation of the error into 
subcomponents). The new result is -3.2 • 10" 10 below 
that of our previous evaluation [10]. This shift is com- 
posed of —0.7 from the inclusion of the new, large photon 
angle data from KLOE, +0.4 from the use of preliminary 
BABAR data in the e+e" 7r+7r"27r° mode, -2.4 from 
the new high-multiplicity exclusive channels, the reesti- 
mate of the unknown channels, and the new resonance 
treatment, —0.5 from mainly the four-loop term in the 
QCD prediction of the hadronic cross section that con- 
tributes with a negative sign, as well as smaller other 
differences. The total error on a h ad TO j g slightly larger 
than that of Ref. [10] owing to a more thorough (and con- 
servative) evaluation of the inter-channel correlations. 

Adding to the result (21) the contributions from higher 
order hadronic loops, —9. 79 ±0.09 [59], hadronic light-by- 
light scattering, 10.5±2.6 [61] (cf. remark in Footnote 8), 
as well as QED, 11 658471.809 ± 0.015 [62] (see also [57] 
and references therein), and electroweak effects, 15.4 ± 
O.lhad ± 0.2 H i gg s [63-65], we obtain the SM prediction 



af 1 = 11 659 180.2 ± 4.2 ± 2.6 ± 0.2 (4.9 t 



(22) 



where the errors account for lowest and higher order 
hadronic, and other contributions, respectively. The re- 
sult (22) deviates from the experimental average, a® xp = 
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FIG. 10: Standard Gfitter electroweak fit result [55] (green 
shaded band) and the result obtained for the new evaluation 
of Aahnd(Mz) (red solid curve). The legend displays the cor- 
responding five-quark contribution, A«[^ d (M§), where the 
top term of -0.72 ■ 10" 4 is excluded. A shift of +7 GeV in 
the central value of the Higgs boson is observed. 

11 659 208.9 ±5.4 ±3.3 [12, 57], by 28.7 ±8.0 (3.6 cr). 8 

A compilation of recent SM predictions for a M com- 
pared with the experimental result is given in Fig. 9. 

Update of T-based g — 2 result. The majority of 
the changes applied in this work, compared to our previ- 
ous one [10], will similarly affect the T-based result from 
Ref. [15], requiring a reevaluation of the corresponding 
r-based hadronic contribution. In the T-based analy- 
sis [67], the tt + tt~ cross section is entirely replaced by the 
average, isospin-transformed, and isospin-breaking cor- 
rected t — > n~ii a v T spectral function, 9 while the four- 
pion cross sections, obtained from linear combinations of 
the t~ — » 7r - 37r°i/ T and t~ — » 2ir~ ir + ir° v T spectral func- 
tions, 10 are only evaluated up to 1.5 GeV with t data. 
Due to the lack of statistical precision, the spectrum is 
completed with e + e~ data between 1.5 and 1.8 GeV. All 
the other channels are taken from e + e~ data. The com- 
plete lowest-order T-based result reads 

a had,LOj T j = ?01 5 ± 3 5 ± Lg ± 2 4 ± q 2 ± 0.3 , (23) 

where the first error is r experimental, the second esti- 
mates the uncertainty in the isospin-breaking corrections, 
the third is e + e~ experimental, and the fourth and fifth 



stand for the narrow resonance and QCD uncertainties, 
respectively. The T-based hadronic contribution deviates 
by 9.1 ± 5.0 (1.8 a) from the e + e~-based one, and the 
full T-based SM prediction a^r] = 11 659 189.4 ± 5.4 
deviates by 19.5 ± 8.3 (2.4 a) from the experimental av- 
erage. The new T-based result is also included in the 
compilation of Fig. 9. 

Running electromagnetic coupling at Af|- The 

sum of all hadronic contributions from Table II gives for 
the e + e _ -based hadronic term in the running of a(M z ) 

Aa had (Mf ) = (275.0 ± 1.0) ■ 10~ 4 , (24) 

which is, contrary to the evaluation of a^ ad ' LO , not dom- 
inated by the uncertainty in the experimental low-energy 
data, but by contributions from all energy regions, where 
both experimental and theoretical errors have similar 
magnitude. 11 The corresponding T-based result reads 
Aa ha d(M|) = (276.1 ± 1.1) • 10~ 4 . As expected, the re- 
sult (24) is smaller than the most recent value from the 
HLMNT group [58] Aa had (Ml) = (275.5 ± 1.4) • 10~ 4 . 
Owing to the use of perturbative QCD between 1.8 and 
3.7 GeV, the precision in Eq. (24) is significantly im- 
proved compared to the HLMNT result, which relies on 
experimental data in that domain. 12 

Adding the three-loop leptonic contribution, 
Aaiep(Aff) = 314.97686 • 10~ 4 [72], with negligible 
uncertainties, one finds 

a- 1 (M 2 Z ) = 128.952 ± 0.014 . (25) 

The running electromagnetic coupling at Mz enters at 
various levels the global SM fit to electroweak precision 
data. It contributes to the radiator functions that modify 
the vector and axial- vector couplings in the partial Z bo- 
son widths to fermions, and also to the SM prediction of 
the W mass and the effective weak mixing angle. Over- 
all, the fit exhibits a —39% correlation between the Higgs 
mass (Mh) and Aahad(Aff) [55], so that the decrease in 
the value (24) and thus in the running electromagnetic 
coupling strength, with respect to earlier evaluations, 
leads to an increase in the most probable value of Mh re- 
turned by the fit. 13 Figure 10 shows the standard Gfitter 
result (green shaded band) [55] , using as hadronic contri- 
bution Aa had (M§) = (276.8 ± 2.2) • 10" 4 [59], together 
with the result obtained by using Eq. (24) (red solid line). 



8 Using the alternative result for the light-by-light scattering con- 
tribution, 11.6 ± 4.0 [66], the error in the SM prediction (22) 
increases to 5.8, and the discrepancy with experiment reduces to 
3.2<7. 

9 Using published t — > Tr~n°v T spectral function data from 
ALEPH [68], Belle [69], CLEO [70] and OPAL [71], and using 
the world average branching fraction [51] (2009 PDG edition). 

Similar to Footnote 2, coarse isospin-breaking corrections with 
100% uncertainty are applied to the four-pion spectral functions 
from t decays [16]. 



11 In the global electroweak fit both cts(Mz) and Aa^^Mi) are 
floating parameters (though the latter one is constrained to its 
phenomenological value). It is therefore important to include 
their mutual dependence in the fit. The functional dependence 
of the central value of Aa^^(Ml) on the value of as(M^) ap- 
proximately reads 0.37 ■ 10 -4 x (<x s {M%) - 0.1193)/0.0028. 

12 HLMNT use perturbative QCD for the central value of the con- 
tribution between 1.8 and 3.7 GeV, but assign the experimental 
errors from the BES measurements to it. 

13 The correlation between Mh and Aahad(jWjj) reduces to —17% 
when using the result (24) in the global fit. 
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The fitted Higgs mass shifts from previously 84_23 GeV 
to 91 GeV. The stationary error of the latter value, 
in spite of the improved accuracy in Aahad(M^), is due 
to the logarithmic Mjj dependence of the fit observables. 
The new 95% and 99% upper limits on M H are 163 GeV 
and 193 GeV, respectively. 

VI. CONCLUSIONS 

We have updated the Standard Model predictions of 
the muon anomalous magnetic moment and the running 
electromagnetic coupling constant at M§ by reevaluating 
their virtual hadronic contributions. Mainly the reesti- 
mation of missing higher multiplicity channels, owing to 
new results from BABAR, causes a decrease of this con- 
tribution with respect to earlier calculations, which — on 
one hand — amplifies the discrepancy of the muon g — 2 
measurement with its prediction to 3.6 a for e + e~-based 
analysis, and to 2.4 a for the r-based analysis, while — on 
the other hand — relaxes the tension between the direct 
Higgs searches and the electroweak fit by 7 GeV for the 
Higgs mass. 

A thorough reestimation of inter-channel correlations 
has led to a slight increase in the final error of the 



hadronic contribution to the muon g — 2. A better preci- 
sion is currently constricted by the discrepancy between 
KLOE and the other experiments, in particular BABAR, 
in the dominant 7r + 7r~ mode. This discrepancy is corrob- 
orated when comparing e + e~ and r data in this mode, 
where agreement between BABAR and t data is ob- 
served. 

Support for the KLOE results must come from a cross- 
section measurement involving the ratio of pion-to-muon 
pairs. Moreover, new 7r + 7r~ precision data are soon ex- 
pected from the upgraded VEPP-2000 storage ring at 
BINP-Novosibirsk, Russia, and the improved detectors 
CMD-3 and SND-2000. The future development of this 
field also relies on a more accurate muon g — 2 measure- 
ment, and on progress in the evaluation of the light-by- 
light scattering contribution. 

We are most grateful to Martin Goebel from the Gutter group for 
performing the electroweak fit and producing Fig. 10 of this paper. 
We thank Changzheng Yuan, Andreas Backer and Claus Grupen 
for information on the radiative corrections applied to BES and 
PLUTO data. This work is supported in part by the Talent Team 
Program of CAS (KJCX2-YW-N45) of China. 
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